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ABSTRACT 

We report on the first detection of GeV high-energy gamma-ray emission 
from a young supernova remnant with the Large Area Telescope aboard the 
Fermi Gamma-ray Space Telescope. These observations reveal a source with 
no discernible spatial extension detected at a significance level of 12.2a above 
500 MeV at a location that is consistent with the position of the remnant of the 
supernova explosion that occurred around 1680 in the Cassiopeia constellation - 
Cassiopeia A. The gamma-ray flux and spectral shape of the source are consistent 
with a scenario in which the gamma-ray emission originates from relativistic 
particles accelerated in the shell of this remnant. The total content of cosmic 
rays (electrons and protons) accelerated in Cas A can be estimated as VFcr — 
(1-4) X 10^^ erg thanks to the well-known density in the remnant assuming that 
the observed gamma-ray originates in the SNR shell(s). The magnetic field in 
the radio-emitting plasma can be robustly constrained as > 0.1 mC, providing 
new evidence of the magnetic field amplification at the forward shock and the 
strong field in the shocked ejecta. 

Subject headings: acceleration of particles — ISM: individual(Cassiopeia A) — 
radiation mechanisms: non-thermal 
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Introduction 



Supernova remnants (SNRs) have long been considered as the primary candid ates for the 



origin of Galactic Cosmic Rays (CRs). Specifical l y, diffusive shock acceleratio n ( lBelllll978 



Blandford fc OstrikedllQTSi : IJones fc Ellisonlll99ll : iMalkov fc O'C Drurvl 120011 ) is widely ac- 
cepted as the mechanism by which charged particles can be accelerated to very high energies 
at collisionless shocks driven by supernova explosions. To maintain the energy density of 
the Galactic CRs, the kinetic energy released in supernova explosions has to be efficientl y 
transferred to CRs with a conversion efficiency of ~ 10% ( iGinzburg fc Syrovatskiil Il964l ). 
However, it has not yet been confirmed whether strong shock waves in SNRs are indeed ca- 
pable of transferring this amount of energy into the acceleration of CR ions, and in doing so 
generating a CR energy density comparable to the energy density contained in the expansion 
ram pressure of the supernova. 

Shock energy can be converted to another form of energy: turbulent magnetic fields. 
Recent theoretical work indicates that the magnetic field can be largely amplified by stream- 



(Bell & Lucek 


2001; 


Bell 


2004) 



have shown that the magnetic field at supernova shocks can be amplified way beyond the 
factor of ~ 4 that is expected for standard compression of interstellar magnetic fields in the 



absence of CR modification of the shock's hydrodynamic s tructure (jUchiyama et al.l 12007 



Volk et al.l l2005l : iBamba et al.l l2005t IVink &: Lamingl |2003[ ). The amplified magnetic field 
would allow acceleration of protons up to the knee in the CR spe ctrum at ~ 10^^ eV, the 



presumed endpoint of the Galactic component (iBell fc Lucekll200ll ) 



Cassiopeia A (Cas A) is the remnant of one of a handful of historical supernovae; the 
explosion around AD 1680 that gave rise to Cas A has probably been the last Galactic 
superno va witnessed by hu mans (lAshworthl Il980t ) . It is the brightest radio source in our 
Galaxy (IBaars et al.l 119771 ) and its overall brightness across the electromagnetic spectrum 
makes it a unique laboratory for studying high-energy phenomena in SNRs. Cas A has an 
angular size of 2.5' in radius corresponding to a physical size of 2.34 pc at a distance of 



3.41q;^ kpc (IReed et al.lll995[). Cas A was the first SNR detected in TeV gamma rays, first 



by HEGRA flAharonian et al.ll200lh and later confirmed by MAGIC flAlbert et al.ll2007l ) and 
VERITAS (jHumenskyll2008l ) . firmly establishing the existence of multi-TeV particles. These 
higher-energy gamma rays are generally attributed to particles accelerated in the shock 
waves of the SNR, even though the TeV gamma-ray source is not resolved. The emission 
mechanism(s) responsible for the TeV gamma-ray emission remained unset tled. EGRET 



repor ted only upper limits on the gamma-ray emission in the GeV band (lEsposito et al. 
1996h . 
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The advent of the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space 
Telescope for the first time makes it possible to detect GeV gamma rays from the shell of 
Gas A. The relativistic bremsstr ahlung of accelerat ed electrons has been predicted to lie 
above the sensitivity of the LAT (lAtoyan et al.ll2000l ). The prediction of the bremsstrahlung 
flux is controlled by the strength of the magnetic field and therefore th e gamma-ray flux 
provi des a measure of the amplified magnetic field in the supernova shell (IGowsik fc Sarkar 



1980l ). Inverse Gompton (IG) scattering is another important mechanism of gamma-ray 



production by high-energy electrons; seed photons that can be upscattered to gamma-rays 
include the interstellar radiation field, the cosmic microwave background (GMB), and the 
far-infrared (FIR) radiation by Gas A itself. In addition to these leptonic emissions, the 
decay of vr^-mesons produced in inelastic collisions between high-energy protons (and nuclei) 
and background gas may contribute to the GeV gamma-ray flux. The importance of the 
GeV observations of young SNRs is emphasized by the fact that the energy density of the 
accelerated particles is measurable thanks to the well-constrained gas and radiation density 
in Gas A in addition to that of the magnetic field. In this paper, we report the discovery of 
GeV gamma-ray emission coincident with Gas A based on the first year observations with 
the Fermi LAT. 



2. Observations 



The LAT onboard Fermi is a pair-conversion gamma-ray detector operating between 
20 MeV and 300 GeV. The LAT has a wide field of view of ~ 2.4 sr at 1 GeV, and observes 
the entire sky every 2 orbits (~ 3 hr f or Fermi'' s orbi t at an altitude of ~ 565 km). The full 



details of the instruments are given in lAtwood et al.l (120091 ). 



A total exposure of 3 x 10 cm s (at 1 GeV) has been obtained for Gas A during 
the period between August 4th 2008 and September 4th 2009, corresponding to 394 days of 
observations. The data analysis was performed usin g the LAT Science Tools package with the 
P6_V3 post-launch instrument response function ( IRando et al.l 120091 ). The standard event 
selection for source analysis, resulting in the strongest background- rejection power {diffuse 
event class) was applied. In addition photons coming from zenith angles larger than 105° 
were rejected to reduce the background from gamma rays produced in the atmosphere of the 
Earth. The analysis was further restricted to the energy range above 200 MeV where the 
uncertainties in the effective area become smaller than 10%. 
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3. Analysis and Results 

The analysis of the gamma-ray event data was performed using gtlike, which is avail- 
able as part of Fermi Science Toolfl The tool gtlike employs a maximum likelihood 
techni que to assess the s tatistical significance of sources and to estimate spectral param- 



eters ( iMattox et al.lll996l ). The background gamma-ray model includes background sources 
from the 11-month catalog of Fermi-LAT sources at fixed coordinates, the galactic diffuse 
emission (using gllJem_v02.fit) and an isotropic component (using isotropic_iem_v02.txt). 
The maximum likelihood analysis is performed inside a region-of-interest (ROI) of 10° ra- 
dius centered on Cas A. The independent tool sourcelike was employed to determine the 
source position and constrain the angular extent of the source. Sourcelike performs a max- 
imum likelihood fit to simultaneously optimize position and extent of the source (given an 
assumption about the shape of the source) in independent energy bands. Sourcelike can also 
be used to assess the test-statistic (TS) value and to compute the spectra of both extended 
and point- sources. The position determined with sourcelike is consistent with what was 
found by the automated LAT catalog tool. The spectral parameters between sourcelike and 
gtlike for this best-fit position are found to be consistent, yielding an independent check of 
the main results of the paper. 

The analysis clearly shows a source above the background coincident with the SNR 
Cas A (see Fig. [T]). The source is detected at a significance level of 12.2 cr (or a TS value of 
148) above the background at a best fit position of 0:2000 = 23h23ml7.5s, ^2000 = 58h49m43.1s 
with a statistical uncertainty on the location of 1.0' (68% confidence level). The systematic 
error in the position due to alignment of the telescope system and inaccurate description 
of the point-spread function of the instrument is estimated to be 0.3'. The position is 
spatially c oincident with Cas A and with the M AGIC and VERlTAS-detected TeV gamma- 



ray source (lAlbert et al.ll2007l : lHumenskyll2008l ) as shown in Fig. [2l An upper limit on the size 



of the gamma-ray emission has been obtained in a maximum likelihood fit by investigating 
the decrease of the likelihood with increasing source size. Under the assumption of a Gaussian 
shape this 68% upper limit amounts to 3.5' (l-a radius). This value is larger than the radius 
of the SNR (2.5'), which thus currently cannot be resolved by the Fermi-LAT . A check for 
variability in this data set with one-month time bins shows no sign for any change in the 
gamma-ray flux with time, indicating a steady source of emission. 

A spectral analysis of the gamma-ray emission has been performed in an energy range 
of 200 MeV-50 GeV. The source is detected only above 500 MeV. The spectrum shown in 



^Software and documentation of the Fermi Science Tools are distributed by the Fermi Science Support 



Center at http://fermi.gsfc.nasa.gov/ssc 
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Fig. |3] reveals a relatively fiat energy distribution and can be fitted by a power law with 
spectral index of F = 2.0 ±0.1 between 0.5 and 50 GeV. The systematic uncertainty in 
the spectral index determination from the uncertainty in the normalisation of the Galactic 
diffuse emission (conservatively estimated to be 10%) amounts to ~ 0.1. The integral photon 
flux above 500 MeV amounts to (8.7 ± 1.3) x 10~^ photons cm~^ s~^. A lilcelihood-ratio 
test was performed to check the presence of a spectral cutoff using a spectral function of 
dN/dE = KE~^ exp{—E/Ecut)- The presence of the exponential cutoff is not statistically 
signiflcant, given the resultant likelihood ratio of —2\n{Lpi^/Lcutos) = 2.6. We formally 
report the best-fit parameters of the cutoff power-law model to make a comparison with 
LAT pulsars: F = 1.7 ± 0.2 and E^^t = 16 ± 9 GeV. 



4. Discussion 



The detection of gamma-ray emission from the direction of Gas A raises the following 
questions: "Is there a compact source such as a pulsar dominating the emission?" and 
"What is the radiating particle population responsible for the emission?". Two emission 
scenarios seem plausible: emission from the central compact object (GGO) in Gas A or 
emission generated in the SNR. Since the GeV gamma-ray source coinciding with Gas A is a 
point-source in the LAT, these two scenarios cannot be distinguished on positional grounds 
therefore circumstantial evidence has to be considered. 



The point-like central X-ray source (ITananbaunj 119991) is generally tho ught to be the 
left-over of the explosion of the massive progenitor star (jPavlov et al.ll2000l ) and known as 
a GGO given that it is radio quiet, un-pulsed in X-rays, and has an X-ray spectrum de- 
scribed by blackbody with characteristic temperatures of about 0.4 keV without indication 
of a non-thermal component. A recent Chandra obse rvation resulted in a 3a limit on the 
pulsed fraction of 16% for periods larger than 0.68 s ( jPavlov fc Lunall2009[ ). The fact that 
no pulsation has been detected at any waveband for the GGO does not rule out, that the 
neutron star is emitting gamma rays. The Fermi-LAT has detected several neutron stars 
as gam ma-ray-only pulsars, pul sars not previously known from observations in other wave- 
bands JAbdo et al 1120081 . l2009ah . 



We searched for gamma-ray pulsations from the source coincide nt with Gas A usin g 
the full data set. We applied the standard time-differencing technique (lAtwood et al.ll2006l ). 
using a maximum frequency of 64 Hz, and a long time-difference window of ~12 days. 
We found no evidence for pulsations which makes the neutron star scenar io less likely. A 
comparison of the blind search pulsars discovered so far (lAbdo et al.ll2009ai ) and the known 
radio pulsars detected by the LAT, suggests that the blind search is approximately 2-3 times 
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less sensitive than a standard pulsation search using the known timii iK solution. This re sults 
in a 5a limit on the pulsed flux of ~ 2 x 10^^ photons cm~^ s~^ (lAbdo et al.ll2009b( ). In 
addition, there is no GeV gamma ray source in the 1st Fermi-LAT source catalog that is 
associated with a known CCO. 

Furthermore, the spectrum does not support a pulsar hypothesis. The energy spectra 
of pulsars are usually flat at energies be low 1 GeV and sh ow exponential cutoffs in the 
energy range between 1 GeV and 8 GeV (lAbdo et al.ll2009bl ). These characteristics do not 
mirror the LAT spectrum which is best described by a power-law shape with an index of 
r = 2.0 ±0.1 up to 50 GeV with no signiflcant sign of a high-energy cutoff. A formal flt with 
an exponential cutoff model yields -Ecut = 16 ± 9 GeV. This is rather high for a gamma-ray 
pulsar; no LAT pulsars show i^cut > 8 GeV so far (lAbdo et al.ll2009a] ). 



The scenario in which the GeV gamma rays are emitted in the shell of the SNR is 
therefore favored. The gamma-ray emission could be produced by electrons accelerated at the 
forward shock through relativistic bremsstrahlung or IC. Alternatively, the GeV gamma-ray 
emission could be predominantly produced by accelerated hadrons through interaction with 
the background gas and subsequent 7r°-decay. Recent studies showed electron acceleration 
to multi-TeV energies is likely to take place also at t he reverse shock in the supernova ejecta 
dUchivama fc Aharonianlboosi Iflelder fc VinMboOsi ). 



First, we consider the shocked circurastella r region between the contact discontinuity 
and the forward shock (see iGotthelf et al.ll200ll ). assumed to have constant magnetic fleld 



of B (a fr ee parameter) and shocked circumstellar medium with a constant density of nn = 
10 cm~^ ( iLaming fc Hwang] l2003l ). Electrons are accelerate d to multi-TeV ene rgies at the 
forward shock as traced by synchrotron X-ray outer fllaments (IHughes et al.ll2000[ ). We adopt 



here an electron acceleration s pectrum Q JE) oc E "^'^^ exp{—E / E^) to match the radio 



IR spectral index of a = 0.67 ( iRho et al.l 120031 ). since both the GeV gamma-ray emission 
and the radio synchrotron emission sample similar electron energies. About half of the 
total radio flux is attributable to this region (so-called plateau), while another half is to 
the reverse shock region (so-called bright ring; see below). Given the radio flux and the 
effective density (riefr = SnjZj(Zj + 1) ~ 26 cm~^), the flux of bremsstrahlung is controlled 
by B. In Figure |31 we show a leptonic model with B = 0.12 mG (red curves), which can 
broadly explain the ob served GeV flux. The maximum energy is set to be E^ = 40 TeV 
( IVink fc Laming] 120031 ). Shown are contributions from bremsstrahlung (dashed) and from 
IC scattering (dotted) produced by accelerated electrons that suffer synchrotron cooling at 
high energies. The bremsstra hlung spectrum cons isting of electron-ion and electron-electron 
components is computed as in lBaring et al.l ( 119991 ). The radiation fleld for the IC component 
is dominated by FIR emission from the Gas A ejecta, characterized by a temperature of 100 K 
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and an energy density of ~ 2 eV cm~^ (IMezger et al.lll986l ) , a factor of 8 larger than the energy 
density in the CMB. The IC/FIR emission exceeds IC/CMB by a factor of 2.7 at a gamma- 
ray energy of 10 GeV. The val ue of -B = 0.12 mG is consistent with B = 0.08-0.16 mG at 



the forward shock estimated by lVink fc Lamingl (120031 ) based on the width of a synchrotron 
X- ray filament. Note h owever that a somewhat higher value of -B ~ 0.3 mG was obtained 
by lParizot et al.l ( l2006l ) using the same filament width but including projection effects. The 
total amount of electrons in this case is We{> 10 MeV) ~ 1 x 10'*^ erg. Also shown in Fig. [3] 
is the case of -B = 0.3 mG (blue curves), which predicts a lower gamma-ray flux than the 
observed one. 

In a scenario in which the gamma-rays are generated by vr^-decay of accelerated hadrons, 
the gamma-ray spectrum can be well matched with either the proton acceleration spec- 
trum Qp{p) oc (a red curve in Fig. H]), or a harder proton spectrum of Qp{p) oc p^"^'^ 
with an exponential cutoff at 10 TeV that is arbitrary introduced (blue curve). Here p de- 
notes momeri t um o f accelerated protons. The gamma-ray spectrum is calculat ed foUowiri g 
Kamae et al.l (120061 ) with a scaling factor of 1.85 for helium and heavy nuclei (IMoril |2009| ). 
The total proton content amounts to Wp{> 10 MeVc~^) ~ 3.8 x 10^^ erg in the case of the 
softer spectrum and to Wp{> 10 MeVc~^) ~ 3.2 x 10^^ erg in the case of the harder proton 
spectrum with the cutoff. In both cases the energy content corresponds to less than 2% of the 
estimated explosion kinetic energy of Egn = 2 x 10^^ erg. Therefore, the cosmic-ray pressure 
would not be large enough to change the hydrodynamics of Gas A. T his is consistent wit h 



Chandra X-ray measurements of the remnant's spatial structure (e.g.. iGotthelf et al.ll200ll ): 
the ratio of the radii of forward and reverse shocks can be reproduced by h ydrodynamical 
models that do not include cosmic-ray acceleration ( iLaming fc HwangI 120031 ). 



Comparing the leptonic and hadronic models, it seems clear that the hadronic scenario 
can better fit the data due to the turnover at low energies that is not well reproduced in 
the leptonic scenario. Given the uncertainties in the diffuse model that have stronger effects 
at the low-energy end, we refrain, however, from strong claims about the radiating particle 
population at this point. Changing the diffuse model normalisation by ±10% (a conservative 
assumption on the uncertainty) largely affects the energy points at and below 1 GeV. The 
resulting effect on the flux point at 1 GeV is a 25% upward and 65% downward shift. A 
more detailed investigation of the lower energy end of the LAT spectrum will be possible 
with future LAT data. 

The shocked ejecta gas that emits strong radio and infrared synchrotron light (known 
as the bright ring) is another potential gamma -ray emitting region. Using Mejecta = 2Mq 
( iWillingale et al.ll2002t iLaming &: HwangI 120031 ) comprised of only oxygen, we obtain Ues = 
noZo{Zo -|- 1) ~ 32 cm~^. This happens to be similar to the value of Ues in the forward 
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shock region. Also, the baryon density of uqAo — 7 cm~^ is close to that in the forward 
shock region. Therefore, the 7r°-decay predictions are essentially the same as those described 
above, though the energy budget is tighter in the reverse shock. The total thermal energy 
stored in the shocked ejecta would be only ~ 1 x 10^° erg. On the other hand, a combination 
of rieff — 32 cm~^ and the large magnetic field in the shocked ejecta {B > 0.5 mG) inferred by 



X-ray variability ( lUchiyama fc Aharonianll2008l ) makes it difficult to attribute the gamma- 



ray emission to the relativistic bremsstrahlung in the shocked ejecta. 

Regardless of the origin(s) of the observed gamma rays, the total content of CRs acceler- 
ated in Gas A can be obtained as Wcr = We + Wp ~ (1-4) x 10^^ erg, and the magnetic field 
amplified at the shock and the field in the shocked ejecta can be constrained as i? > 0.12 mG. 
Even though Gas A is considered to have entered the Sedov phase, the total amount of GRs 
accelerated in the remnant constitutes only a minor fraction (< 2%) of the total kinetic 
energy of the supernova. The bremsstrahlung spectrum and 7r°-decay spectrum have rather 
different predictions below 1 GeV. The hard spectrum below 1 GeV would favor the vr'^-decay 
origin, though the current LAT data quality does not rule out the bremsstrahlung model. 
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France, ASI and INFN in Italy, MEXT, KEK, and JAXA in Japan, and the K. A. Wallen- 
berg Foundation, the Swedish Research Gouncil and the National Space Board in Sweden. 
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Fig. 1. — Smoothed gamma-ray count map of the region surrounding Cas A with side-length 
12°, binned in square pixels of side-length 0.05°. Smoothing was done with a Gaussian of 
width 0.1°. Overlaid are test statistics contours (TS = — 21n(— L)) for the assumption of a 
point-source located at the position of the pixel. Contours of 25, 100, 200, 300, 400, 500, 
600 and 700 TS levels are shown. Only events with energies larger than 800 MeV were used. 
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Fig. 2. — The VLA 20 cm radio map of the Cas A supernova shell (jAnderson fc Rudnick 



19951 ). Shown is the region corresponding to the green square in Figure [TJ The graph shows 



the GeV position with error bars (conservatively adding statistical and systematic errors in 
quadrature) as a cyan circle and the position of the CCO as a y ellow star. Also s hown in the 
plot are the best-fit positions for MAGIC jAlbert et all bon?! ) and VERITAS (jHumensky 



20081). 
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Fig. 3. — Energy spectrum of Cas A in a leptonic emission model. Shown is the Fermi-LAT 
detected emission ( filled circles) in co mparison to the energy spectra det ected by MAGIC 
(black open circles; lAlbert et al.l 120071 ) and VERITAS (green open circles; iHumenskyI |2008| ) . 
The red curves show a leptonic model calculated for B = 0.12 mG while the blue curves show 
one for B = 0.3 mG. Dashed lines show the bremsstrahlung components with Ues — 26 cm~^, 
and dotted lines show the IC component. 
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Fig. 4. — Same as Fig. |3]but in a hadronic emission model. Shown are 7r°-decay spectra for 
two possible proton spectra. The blue line shows a harder proton index of 2.1 with a cutoff 
energy of 10 TeV {Wp = 3.2 x 10^^ erg for = 10 cm~^). The red line shows a softer proton 
index of 2.3 without a cutoff. 



